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Abstract 

We present high statistics simulations for 2-d percolation clusters in the "bus bar" 
geometry at the critical point, for site and for bond percolation. We measured their 
backbone sizes and electrical conductivities. For all sets of measurements we find large 
corrections to scaling, most of which do not seem to be described by single powers. 
Using single power terms for the corrections to scaling of the backbone masses, we 
would obtain fractal dimensions which are different for site and bond percolation, 
while the correct result is = 1.6432 ± 0.0008 for both. For the conductivity, the 
corrections to scaling are strongly non-monotonic for bond percolation. The exponent 
t' = t/v is measured as 0.9826 ± 0.0008, in disagreement with the Alexander-Orbach 
and other conjectures. 

1 Introduction 

The critical behavior of percolation is quite well understood, mainly because of its relation- 
ship to the one-state Potts model Jl|. Due to this relationship and conformal invariance, 
all "thermal" exponents are exactly known in d = 2. There are however a number of crit- 
ical exponents for percolation which have no thermal analogue. The most important of 
these are the backbone dimension and the conductivity exponent t which is defined by 
° ~ (p ~ PcY above the critical point. Here and in the following we assume a 'bus bar' 
geometry of size L x L, and a is the global conductivity, with each empty site/bond being 
an insulator and each occupied site/bond having a constant finite resistivity. In a bus bar 
geometry two opposite sides of a quadratic lattice are connected to superconductors, while 
the other two sides are open. The backbone is defined as the set of lattice sites/bonds which 
are connected to both bus bars through mutually non-intersecting paths. Related to Db 
and t are various exponents for conduction exactly at the critical point, for conductivity 
of conductor/superconductor mixtures, and for random walks starting either on the infinite 
incipient cluster, on the backbone, or at arbitrary lattice sites ||. Instead of t we shall in 
the following discuss mostly t' = tv, with v = 4/3 being the correlation length exponent. 



Exactly at threshold, the conductivity of a lattice of size L scales as a ~ L 1 . At present, 
the most precise Monte Carlo simulations give t' = 0.9745±0.0015 and D b = 1.647±0.004 
|3j. These seem to be more precise than estimates based on exact enumerations || or other 
numerical methods. 

Theoretical attempts to compute these "athermal" exponents have turned out to be less 
easy. There have been numerous conjectures, most of them rather ad hoc. The best known is 
due to Alexander and Orbach ||, and gives in 2 dimensions t' = 91/96 = 0.94791 .... It was 
for some time considered as very accurate, but if we accept the above numerical estimate, it 
is ruled out by about 18 standard deviations. 

It seems that there are no conjectures for t' based on conformal invariance, but there are 
several conjectures for D\,. In two independent attempts, Larsson [[/] and Saleur || obtained 
the value Dj, = 25/16 = 1.5625 which is obviously in serious disagreement with the above 
numerical estimate (and with other recent numerical estimates || [TO])- More recently, Huber 
obtained Df, = 79/48 = 1.64587. This is in excellent agreement with the simulations, 
but its theoretical justification is not very clear. 

The present investigation was started for several reasons. The first is that the random 
number generator used in || (the Stoll-Kirkpatrick generator R250 JT2|) had obviously cre- 
ated some problems. For L — ► oo the spanning probability should tend to 1/2 in the square 
bus bar geometry, with corrections of size 1/L |TB|. While the convergence to 1/2 was 
seen in 0, it was much slower than it should have been. As a consequence, the value of 
quoted above should be taken with some caution. Similar problems with R250 were found 
in other depth-first algorithms |13], |14|| . The spreading simulations presented in sec. 2 of |4j] 
should be unaffected by these problems since they used a breadth-first algorithm for which 
R250 seems to be safe. 

Secondly, we considered the situation concerning the conductivity exponent as unsatis- 
factory. The random number generators used in |3|] were obviously even worse than R250. 
The data for site percolation presented in table 1 of show fluctuations which are much 
larger than the supposed statistical errors, and could hardly be blamed on corrections to 
scaling, because of their irregular behavior. In addition, the estimates obtained in |3[ for site 
and bond percolation are not compatible within their error bars. 

Finally, it seemed that many of the papers in the recent literature used quite slow al- 
gorithms, while fast algorithms have been applied only with moderate statistics. This is 
particularly true for conductivity. There, the fastest algorithm by far (for the square lat- 
tice!) seems to be the one by Lobb and Frank [|17], |j~8"| . Its time complexity is roughly 
L 2 logL. This should be compared to the time complexity L 3 x N for strips of size L x N 
(with N ^> L) for the algorithm used in || . Nevertheless, while several months of CPU time 
on a special purpose computer were used in ||, only a few hours on an IBM 3081 mainframe 
were spent on the published data for the Lobb-Frank algorithm. [] In addition, there exists 
a vast literature in which random walks were used to estimate t via the Einstein relation |2j . 
I implemented several variants (walks on percolating clusters only, walks on backbones only, 
walks on all sites). Even though I spent substantial amounts of CPU time on them, results 
were extremely poor compared to the results described below. 

1 Dr. Frank has, however, produced a substantial amount of yet unpublished data fl9| . I am greatly 
indebted to him for sending me these data. They were extremely helpful in convincing me that my own 
simulations were correct. 



Similarly, although the algorithm used for backbone identification in Q has a complexity 
L D f where Df = 91/48 = 1.895..., much slower algorithms are still in use. This seems 
particularly true for various versions of the burning algorithm [EH], I22J, 123, 1231, 0, 551. An 



algorithm with complexity slightly larger than L 2 was used recently in |H| |20|. It should be 
pointed out that the algorithm of [f§] works only for strictly planar graphs (no such restriction 
holds for burning and for the algorithms of [10, 20|), but a slightly more complex algorithm 
with essentially the same asymptotic behavior was given long ago by Tarjan for arbitrary 
graphs (Tarjan's algorithm finds the 'biconnected' part of any graph which for percolation 
is just the backbone). 

Indeed, the fastest algorithm for estimating conductivity turned out to be a combination 
of the Lobb- Frank algorithm with the backbone-identification algorithm of 0]. First the 
backbone is found (time L Df ), and then the Lobb- Frank algorithm is run on the backbone. 
The time for the latter is still > 0(L 2 ), not L Db logL, as one might have expected naively, 
and one has to add the CPU time needed for finding the backbone. Nevertheless, total 
absolute CPU times are reduced by roughly a factor 4 for the largest values of L studied 
here. 

In retrospect, however, the main result of the paper is something unexpected, namely the 
very strange behavior of the corrections to scaling. Field theory predicts that the leading cor- 
rections are power behaved, with universal exponents derivable by means of renormalization 
group techniques. As a consequence it is often assumed that a single power behaved correc- 
tion term is sufficient to describe all deviations from the asymptotic scaling laws. In many 
cases this is correct, in particular if the data which are to be interpreted are not too precise. 
But there are counter examples. The best known maybe is the spherical model where Luck 
|27j showed analytically that corrections to scaling are non-monotonic. Another example is 
spreading of 3-d percolation where Monte Carlo simulations would give non-universal critical 
exponents if these simulations were fitted by simple power-behaved corrections p8| . 

The corrections to scaling discussed below seem even more strange, in particular those 
for the conductivity. A priori one should expect that bond percolation is better behaved 
than site percolation |fl7| , |j~8| . The first reason is that the critical point is exactly known to 
be p c = 1/2 |1J. Secondly, if the geometry is properly chosen (together with the bus bars, 
the lattice has shape L X (L + 1) [0), duality shows that the spanning probability (i.e. the 
probability that the backbone is non-empty and that o > 0) is exactly 1/2 at p = p c , for all 
L. Indeed, corrections to scaling are smaller for bond percolation than for site percolation. 
But while they are reasonably well described by a single power for the latter, they show at 
least one full oscillation for bond percolation. Indeed, from the bond percolation data alone 
we cannot exclude the possibility that they continue to oscillate with log L (such log-periodic 
oscillations do occur in a similar conduction problem for large applied voltage [[^J). It is 
only the comparison with site percolation for very large L which makes ongoing oscillations 
unlikely. 

In the next section we will discuss backbones, in sec. 3 conductivities. Our conclusions 
will be drawn in sec. 4. 



2 Backbones 



The backbone of a (site) percolation cluster which spans between two sets A and B of sites 
is defined as the set of those sites which are connected to both A and B by two mutually 
non-intersecting paths. Except for 'Wheatstone bridge' type configurations, these are also 
those sites which carry non-zero current if the lattice is made of a finite resistance conductor, 
and a potential difference is applied between A and B. In the geometry used in the present 
paper, A and B are the bus bars which are supposed to have infinite conductivity. 

The backbone is essentially (if an additional bond is added between A and B) what 
is called the set of biconnected nodes in computer science. For general graphs, the fastest 
known algorithm to generate it is a recursive depth-first algorithm [25| which runs in time 



O(N) for a graph of N nodes, provided the number of edges meeting at each node is finite. 
Nearly the same iV-dependence was reached recently in a novel algorithm [[Kj which can 
also be extended to rigidity backbones |2(J . At least for large N this is much faster than the 
'burning' type algorithms used traditionally by physicists [[H], |22|, ^3], |4], |9|, |25| . For strictly 
planar graphs, an even faster algorithm was given in M and applied to 2-d percolation in the 
'bus bar' geometry (this algorithm cannot be applied, e.g., to 2-d percolation with periodic 
lateral boundary conditions, since this is not a strictly planar graph). This algorithm has 
the same asymptotic complexity as Tarjan's algorithm, but it is roughly twice as fast and 
uses about half of the memory, since it needs one data structure less and needs only one pass 
through all sites, instead of two passes in Tarjan's algorithm. 

For site percolation we used essentially the same algorithm as in || , where also a detailed 
description can be found. The main difference is that we realized that we do not have to 
erase disorder configurations (i.e., set all sites variables again to some default value) before 
starting a new configuration. Instead, we attach a long integer Sj to each site, and initialize 
all Si to zero before starting the simulations. While building the first configuration, we 
consider a site as untested if Sj < sf~\ and increase Sj after testing it in the same way as 
in [|J, except that the information whether the actual path is upward or downward is no 
longer encoded in the sign of s^ (as in M), but in an additional array. After completing 
the fc-th configuration, s- fc+1 ^ was set equal to maxj Sj, and the process was repeated. For 
bond percolation the algorithm is slightly more complicated, but is a rather straight forward 
modification. For both cases, the algorithms were carefully checked by visually inspecting 
a large number of examples. Another test consisted in comparing the conductivity of the 
backbone with that of the entire lattice. The two should be identical, and were so in all 
tested cases. 

For site percolation we place the bus bars at the x-axis and at the line y — L + 1, and 
use as 'wettable' lattice sites those with < x, y < L. If the backbones defined in this way 
contained N sites, their dimension is defined as 

(N) ~ L Db . (1) 

For bond percolation we place the bus bars at y = and at y = L. Horizontal wettable 
bonds are those with 1 < x < L, < y < L, and vertical wettable bonds are those with 
1 < x < L, < y < L. The backbone can now be defined again as the set of biconnected 
sites, or as the set of biconnected bonds. We call the cardinalities of these sets N s and Nf,. 
For both of them we expect the scaling law eq. ([!]), with the same value of A,. 



As random number generators we used the multiplicative rule i n+ \ = i n *13 13 +l (mod 2 63 ) 
and the four-tap shift-register generator of Ziff |3(J with period 2 9689 — 1. Both gave results 



in perfect agreement. For bond percolation at p = 1/2 the spanning probability was 1/2 
for all values of L. For site percolation we used the value p c = 0.592746 found in fl5 |. The 



spanning probability now converged to 1/2 with a power close to —1, in perfect agreement 
with |Tj| . In order to check for systematic errors due to a wrong estimate of p c , we made 
also runs at p — 0.59276 which is about 30 standard errors larger than the central value of 



15]. These runs showed that the estimates given below should still be correct even if the 



errors in p c were underestimated in [15| by a factor of 2. 



For site percolation, raw data including sample sizes and CPU times are given in table 1. 
Lattice sizes studied were powers of 2, with L up to 4096. In addition we simulated lattices 
with L = 3, 6, 12, . . . 96. For each L < 48 the number of configurations was > 10 8 . This 
number decreased for larger L, see table 1. All simulations were done on Sun Ultra and DEC 
Alpha work stations, but for convenience CPU times are quoted only for the DEC Aplha 
machines. Determining the backbone (and the entire spanning cluster) on a 4096 x 4096 
took in average less than 3 seconds. For larger L the CPU time increased slightly faster 
(roughly ~ L 2 ), probably due to increasingly frequent cache misses. 

For bond percolation we used similar lattice sizes and went up to similar statistics. 
Details are given in table 2. CPU times were not measured separately, as all these runs were 
made together with conductivity measurements, but they should be similar to those for site 
percolation. 

In order to obtain estimates for Db, we computed effective dimensions according to 

]og[N(2L)/N(L/2)] 
6 ' cff( } = fo£i • { ) 

Here N stands for the number of sites in site percolation, for the number of sites in bond 
percolation, or for the number of bonds in bond percolation. This gives us three sequences 
of effective dimensions. 

To obtain a first rough impression, we plot them in fig.l against 1/L. Only values for 
L > 20 are shown. All three curves seem to converge to the same value ~ 1.643, but with 
very different corrections to scaling. 

For site percolation (lowest curve) we see the strongest corrections. Here -D 6 eff increases 



strongly with L, explaining why was underestimated in the first analyses |2T| , P2]] . The 
data seem to follow a smooth line, suggesting a single dominant correction term. Plotting 
the data against 1 / L A and looking for a straight line does not produce entirely satisfactory 
results (showing that more than one term is needed), but if one insists in a single term, its 
power is A « 0.67. The quality of this fit can be judged from fig. 2 which shows the effective 
exponent after subtracting a term a/L A with A = 0.67 and a such that the data are fitted for 
L > 32. Notice that we show a larger range of L in fig. 2, and the fit with a single correction 
term becomes indeed bad for L < 24. 

The same kind of analysis was also done for bond percolation. Again we see that the 
data plotted against 1/L A do not give perfect lines for any A, but the best fitting values 
are definitely larger than for site percolation: A = 1.11 for N s , and A = 1.30 for iV 6 . As 
expected, it is bond percolation with Nj, defining the backbone mass which gives the largest 
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Figure 1: Effective backbone dimensions, estimated according to eq.(^), plotted against 1/L. 
The lowest curve is for site percolation, the uppermost for sites belonging to bond percolation 
backbones, and the middle for bonds in bond percolation backbones. When error bars are not 
visible, they are smaller than the sizes of the symbols. 

A and, as seen from fig.l, also the smallest amplitude for the correction to scaling. The 
results obtained by subtracting the fitted correction terms are also plotted in fig.2. 

If the above procedure were physically meaningful, the three curves in fig.2 should be 
horizontal and should collapse. While the former is at least true for two of them (the curve for 
backbones denned via bonds in bond percolation, which should show cleanest behavior, is not 
flat), the latter is definitely not true. The subtraction reduced considerably the differences 
between the curves. But the error bars are so small that the remaining differences still are 
highly significant. 

Indeed it is not surprising that the above procedure does not give satisfactory results. 
From universality one should have expected that A is the same for all three observables. Our 
findings that A varies between 0.67 and 1.20 shows that these can only be effective exponents, 
obtained by fitting by a single power a function which actually is a superposition of several 
powers. Unfortunately, trying to fit the data with two terms, -L\ e ff(L) = Db + a/L Al + b/L A2 
gives either no satisfactory fits (if too many values of L are used in the fit) or ambiguous fits 
(if only very large values of L are used). 

Our final procedure was the following. First of all, the bond percolation data strongly 
suggest that the leading non-analytic corrections have exponent > 1. The leading term is 
then analytic, ~ 1/L, and the third important term is again analytic, ~ 1/L 2 . We fitted the 
data therefore with ansatzes 



L\ eff (L) = D b + a/L + b/L A + c/L 2 



(3) 



with the same exponent A for all three data sets. The fit (using all L > 10) gave A = 1.42. 
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Figure 2: Effective backbone dimensions, after subtracting terms a/L A , with a and A obtained 
by fitting -D^eg to ansatzes -D^eff = -D& + a/L A (these fits are done separately for each of the three 
curves of fig.l). Notice that the range on the y-axis has been reduced substantially, but the three 
curves do not collapse. In particular, they do not extrapolate to the same value for L — > oo. The 
dashed horizontal line is the prediction of [11 1. 



We do not quote any error bars since this value might change considerably if we include 
even more terms, or if we change the range of L used for the fit. Also the constants a, b, 
and c should not be meaningful. But we checked that the value of D b obtained in this way 
is very robust and agrees for all three data sets, see fig. 3. Figure 3 shows some deviations 
from horizontal lines for very large L, but they seem to be statistical fluctuations. Our final 
estimate is 

D b = 1.6432 ± 0.0008 . (4) 

The error bars are subjective. They were estimated by comparing with similar figures ob- 
tained by changing the range of L and including a fourth term with A = 2.5, but restricting 
its amplitude to a value < |c|. 

Our estimate is five times more precise than that of |3| , but agrees with it within the error 
bars. It is of course completely incompatible with the prediction of |7|, [8|, but it disagrees 
also with Huber's prediction by about 3 standard deviations. It is just outside the error bars 
of the estimate 1.650 ± 0.005 of [|IJ. 

3 Conductivity 

For conductivity measurements we used the Lobb-Frank algorithm, after removing all sites 
resp. bonds not belonging to the backbone. Since implementation of this algorithm is not 
quite trivial, and since our results were rather surprising, we performed very extensive tests. 
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Figure 3: Similar to fig. 2, but with three correction terms subtracted from each curve. Two of 
these terms are analytic (with exponents —1 and —2), the third has a common exponent A = —1.42. 



We compared the results of the Lobb-Frank algorithm with several other algorithms. In 
particular we used several variants of Fogelholm's method |5T] and the direct solution of the 
Kirchhoff equations. In all cases we compared also the conductivity of the backbone with 
that of the entire lattice, finding perfect agreement within the numerical accuracy. During 
these checks we also verified that the Lobb-Frank algorithm is by far the fastest. 

In the Lobb-Frank algorithm one imagines wires to be attached to the lower left and 
upper right corners of the lattice. One then uses alternatingly the star-triangle and the 
triangle-star relations to "push" the upper and left boundary of the lattice downward and 
to the right in such a way that the total conductivity is unchanged. When the calculation 
is finished, the whole lattice consists of only the lower and right boundaries. The lower 
boundary will be superconducting, and the entire resistance comes from the right boundary 
which consists of a string of resistors connected in series. 

When removing parts that do not belong to the backbone, one still has to push the 
boundary through the entire lattice, whence the complexity cannot be smaller than 0(L 2 ). 
But most of the star-triangle and triangle-star relations will be trivial, leading to the sub- 
stantial numerical improvement mentioned above. For bond percolation on a lattice with 
L = 4096, the average CPU time on an DEC Alpha with 433 MHz was ca 16 sec per lattice. 
This included the determination of the backbone and, if this was non-zero, the computation 
of the conductivity. 

For bond percolation we computed the conductivity for all configurations. For site per- 
colation we did it only for about half of the configurations, except that we did not compute 
conductivities for site percolation lattices with L = 4096. In all cases we produced data only 
at the exact (resp. numerically precise) percolation threshold. 

Altogether our sample involved more than 10 14 wetted sites. This should be compared 
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Figure 4: Conductivities on lattices of size L x L, multiplied by L°.98 and by an arbitrary 
constant. The upper curve is for bond, the lower for site percolation. The bond percolation data 
have been multiplied by a factor 1.2845 in order to coincide for large L with the site percolation 
data. 



to ca. 10 11 wetted sites in 0. In addition, most of our data were obtained for lattices 
with much larger (transverse) size than those of |§ , and the algorithm used in || scales very 
unfavorably with lattice size. It should finally be pointed out that the Lobb-Frank algorithm 
can also be used (with slight modifications) for the strip geometry used in |§, provided the 
periodic lateral boundary conditions are replaced by helical ones. 

Results are shown in fig.4. In this figure we present both site and bond percolation 
conductivities, both multiplied by L 098 . On the x-axis is plotted logL. While we see a 
single-humped curve for site percolation, the data for bond percolation (which should be 
more clean according to the theoretical arguments given in sec.l) show an additional bump 
with maximum at L — 3. For large L both curves seem to converge, suggesting that the 
leading correction to scaling term is the same for bond and site percolation. 

The bump seen for bond percolation at small L is definitely not a statistical fluctuation 
(the error bars for L < 50 are much smaller that the sizes of the points). It is not due 
to a programming error either, as seen by comparing with the data of |17|, and with 
unpublished data by D.J. Frank |L9| (for L < 4 we also computed a by exact enumerations, 
and obtained results in perfect agreement with the simulations). It is completely unexpected, 
and its origin is theoretically not understood. 

The bond percolation data shown in fig.4 could equally well be fitted by a log-periodic 
oscillation, 

a = L°- 98 (a + 6sin(alog(L/L ))) . (5) 

It is mainly the good agreement with site percolation at large L which makes such an 
extrapolation very unlikely. 
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Figure 5: Effective exponents computed from the data shown in fig. 4, plotted against 1/L 
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The best estimate for the exponent t' is obtained by again defining an effective exponent 
t' eS (L) by means of eq.(|2|), and plotting it against a suitable power of 1/L. If there is a 
single dominant correction to scaling term, we obtain a straight line extrapolating to t' . 
Such a straight line cannot of course be expected for the bond percolation data, in view of 
the oscillation seen in fig. 4. But the site percolation data do show a beautifully straight 
line when plotted against 1/L - 8 , see fig. 5, suggesting A = 0.8. This should of course not 
be taken too serious in view of our experience with the backbone dimension. Therefore we 
base our final estimate of t' not only on the straight line extrapolations in fig. 5, but also on 
similar extrapolations for plots with A ranging from 0.7 to 1.3, and on fits with up to three 
correction terms with the same exponents as in the last section. The error bars in our result 



t' = 0.9825 ± 0.0008 



t = t'u = 1.3100 ±0.0011 



(6) 



include the associated uncertainty. 

This estimate is nearly by a factor 2 more precise than the estimate of || , which itself was 
by far the most precise previous estimate. It is barely compatible with it. It is compatible 
with the estimate of [pl|. It excludes the Alexander-Orbach conjecture, the conjecture t' — 1 
of p2| , and the recent prediction t' = 0.995 ± 0.001 of |33[] . When translated into exponents 
for diffusion in disordered lattices ||, it gives exponents which are more precise than all 
previous estimates. For random walks on the infinite incipient cluster it gives e.g. the 
fractal dimension || 
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(7) 



Direct estimates of d w (not using the Einstein relation) have errors which are larger by at 
least one order of magnitude M. 



4 Conclusion 



We have presented high statistics simulations of 2-d percolation, using what we believe to 
be the fastest algorithms known at present. Our efforts were concentrated at "a-thermal" 
exponents which cannot be computed easily by mapping percolation onto the 1-component 
Potts model and using the standard machinery for spin models. We obtained what we 
believe to be the most precise estimates for the backbone dimension and for the conductivity 
exponent t. Indeed, our data sets are much larger than previous ones, but our estimates for 
critical exponents are not so much more precise because of very substantial corrections to 
scaling. 

For most observables we found that these corrections to scaling could not be described by 
a single term. In some cases (backbone dimension) a single term would superficially seem to 
give a decent fit, but such fits would be rather misleading. In that case it was the comparison 
between bond and site percolation which allowed us to reach a more definite and satisfactory 
conclusion. 

For conductivities, the main phenomenon was an unexpected anomaly for bond percola- 
tion which had a priori been expected to yield the cleanest signal. By itself, this anomaly 
could have been interpreted in various ways, including logperiodic oscillations. Again it was 
only the comparison with site percolation which ruled out such exotic possibilities, and which 
allowed us to obtain a precise estimate of the critical exponent. 

Although the precise values of the critical exponents should be interesting by themselves, 
we believe that the more important message of the present paper is that one should be 
extremely careful with corrections to scaling. This is particularly true as data get more and 
more precise. A simple least square fit on a log-log plot might be appropriate for extremely 
crude data, although it is dangerous even then. Much better is of course the widely used 
strategy of making several least square fits, excluding more and more data points which might 
be outside the proper scaling region. An alternative is to fit ansatzes which contain suspected 
correction terms. This is certainly strongly advised, but as seen from the above examples 
it can also be very misleading, if the corrections to scaling have unexpected structures. 
The main point we want to stress is that the latter might be more common than is often 
appreciated. 
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Table 1: Raw data for site percolation. The second column contains the total number of 
configurations, while column nr. 6 gives the number of configurations for which also o 
was computed. Column nr. 5 gives the relative rms. width of the distribution of N, the 
error on N is given by this number divided by the square root of the number of config- 
urations and multiplied by N. CPU times are measured in hours, and are quoted for a 
433 MHz DEC Alpha workstation. The error on the spanning probability P sp an is given by 

\J ' -N -Pspan (1 -fspan ) • 
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Table 2: Raw data for bond percolation. Column 8 gives total CPU times from runs on 
various Sun Ultra and DEC Alpha machines with different speeds. For L < 4, the second 
column shows the total number of distinct configurations, and all results were obtained by 
exact enumerations. For further explanations see table 1. 



